Orbital-selective Mott transitions in two-band Hubbard models 



Nils Bliimer , Carsten Knecht^ , Krunoslav Pozgajcic^' , Peter G. J. van Dongen 

'^Institute of Physics, Johannes Gutenberg University, 55128 Mainz, Germany 
^Institute for Theoretical Physics, J. W. Goethe University, 604 38 Frankfurt/Main, Germany 

Received 7 June 2006 



Abstract 

The anisotropic two-orbital Hubbard model is investigated at low temperatures using high-precision quantum Monte Carlo (QMC) 
simulations within dynamical mean-field theory (DMFT). We demonstrate that two distinct orbital-selective Mott transitions 
(OSMTs) occur for a bandwidth ratio of 2 even without spin-fiip contributions to the Hund exchange, and we quantify numerical 
errors in earlier QMC data which had obscured the second transition. The limit of small inter-orbital coupling is introduced via a 
new generalized Hamiltonian and studied using QMC and Potthoff 's self-energy functional method, yielding insight into the nature 
of the OSMTs and the non-Fermi-liquid OSM phase and opening the possibility for a new quantum-critical point. 

PACS: 71.10.Fd; 71.27.+; 71.30.+h 
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■ 1. Introduction 

Recently, the suggestion [1] that the two consecu- 
tive phase transitions experimentally observed [2] in 
. Ca2-a;Sra;Ru04 should be interpreted as "orbital-selective 
' Mott transitions" (OSMTs) sparked a cascade of related 
. theoretical papers (see [3] and references therein) . 

Microscopic studies of OSMTs usually consider the 2- 

■ band Hubbard model H = Hi + H2 + H3 , where 

' Hi = — y ^ ^m^ima'^jma + ^ ^ ^ nim^niral , 
{ij)ma irn 

,iU'-S,,,J.)n,i,n^,^, 

include hopping between nearest-neighbor sites i,j with 
amplitude tm for orbital m S {1,2}, intra- and inter- 
] orbital Coulomb repulsion parametrized by U and U' , 
< respectively, and Ising-type Hund's exchange coupling; 

^ima = 4ma<^^ma Cr G {T, i}- In addition, 

contains spin-flip and pair-hopping terms (with 1 = 2, f =| 
etc.). As in [3] , we refer to Hi + H2 + H3 with the isotropic 
coupling = = J as the J-model and to the simpliflcd 
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Fig. 1. a) Low-frequency analysis of self-energy. Main panel: 
weight of singularity at a; = for narrow/wide band (filled/open 
symbols) as determined by polynomial fits [shown in inset for 
U = 2.8,2.6,2.4,2.2,2.0, 1.8 (top to bottom)], b) Spectra for metal- 
lic, orbital-selective, and insulating phases. 

Hamiltonian Hi + H2 as the Jz-model; unless noted, U' = 
U/2, J = U/A so that U' + 2J = U. 

Early theoretical studies suggested that the expected 2 
distinct OSMTs occur only in the full J model [4], but not 
in the Ising type Jz-model [5,6]. Thus, it seemed as if spin- 
flip and pair-hopping terms were essential ingredients to 
orbital-selective physics. However, this is not actually the 
case, as shown in a low- frequency analysis of high-precision 
QMC data [3] in Fig. la: a singularity develops at = 
for Uci ~ 2.0 in the narrow-band self-energy, but only at 
Uci ~ 2.5 in the wide-band self-energy (with linear increase 
for U > Uci/2)- Corresponding spectra (Fig. lb) illustrate 
the characteristics of the 3 distinct phases. 
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Fig. 2. a) QMC estimates of quasiparticle weights Z at T = 1/32: 
a) high-precision data [3] (cirles, crosses) clearly shows kinks at 
Uci ~ 2.0; b) suitable analysis reveals a second kink at Uc2 ~ 2.5. 
The second transition is lost in the noise of earlier data [5] (squares) 
with errors exceeding 100% at both transitions (c). 




Fig. 3. Quasiparticle weights for generalized two-orbital model: at 
T = 1/40, the wide-band transition remains first order (with hys- 
teresis) for small enough inter-orbital coupling o < 0.2. 

Discrete estimates Z = [1 - ImJ^iinT) / (ttT)]-'^ for the 
quasiparticle weight clearly show (only) the narrow-band 
transition (Fig. 2a); a second (wide-band) transition is 
visible as kink only in high-precision QMC data (circles, 
crosses) and after adding a Hnear term (Fig. 2b). Possible 
signals from this transitions cannot be distinguished from 
numerical errors in earlier QMC data [5] (squares) since 
they exceed 100% at both transitions (Fig. 2c). 

2. Limit of small inter-orbital coupling 

Clearly, the resolution of two orbital-selective Mott tran- 
sitions with critical interactions differing by only about 20% 
is a very challenging task, in particular at temperatures at- 
tainable using QMC. In this situation, much insight can be 
gained by abandoning the constraint U' + 2J = U fulfilled 
in nearly all earlier studies and instead studying the limit 
of small inter-orbital coupling. Hence, we consider H = 
Hi + aH2 with < a < 1 so that a = corresponds to un- 
coupled orbitals and a = 1 to the case studied previously. 

It is a priori clear, that for a = each orbital should un- 
dergo a usual Mott transition at an interaction determined 
by the corresponding bandwidth [W — 2 for the narrow, 
= 4 for the wide band); note however, that even in this 
case the QMC results for both orbitals have no scaling re- 
lation at fixed T > 0. This is seen in Fig. 3: For a = 0, a 
large hysteresis region appears in the wide-band quasipar- 
ticle weight (uppermost, grey line) while only a single coex- 
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Fig. 4. a) Discrete QMC estimates of wide-band quasiparticle weights 
for a = 0.2. In OSM phase (2.3 <U < 3.9), Z^ide decays to for 
T — > 0. b) In contrast, the self-energy is practically T-independent 
in QMC (squares: T = 1/20, circles: T = 1/40), consistent with DIA 
for T = (solid lines). 

istcncc point is resolved for the narrow band (at T = 1/40) . 
Very importantly, the wide-band transition evidently re- 
mains first order at small, but significant inter-orbital cou- 
pfing [a = 0.1, a = 0.2). It may be expected that the first- 
order range (in a) increases at lower temperatures which 
suggests that the wide-band OSMT might be very weakly 
first order even at a = 1. However, the alternative of a 
quantum phase transition at some critical value a = ac is 
equally interesting and warrants further investigation. 

3. Nature of orbital-selective Mott phase 

As shown above, the (discretely estimated) quasiparticle 
weight Z is not well suited for detecting the second OSMT. 
In fact, it is even misleading in the OSM phase: as seen in 
Fig. 4a, ^widc is nearly constant as a function of T in the 
metallic Fermi-liquid phase {U = 2.0). However, this ob- 
servable decays to for T ^ in the OSM phase. Naively, 
one might conclude that the OSM phase becomes indistin- 
guishable from the insulating phase for T = 0. However, 
Fig. 4b proves that this is not the case: the self-energies 
are practically T-independent both in the metallic and in 
the OSM phase [with agreement between QMC results for 
different T and self-energy functional theory in dynamical 
impurity approximation (DIA)]; thus, the T dependence in 
Zw\dc is an artifact of the discrete approximation. The non- 
Fermi-liquid character of the wide band in the OSM phase 
is clearly seen as the finite limit of Im S for lo — > i0+. 
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